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Abstract 

We  report  a  systematic  measurement  of  the  space  charge  effect  observed  in  the  few-ps  laser 
pulse  regime  in  laser-based  solid-state  photoemission  spectroscopy  experiments.  The  broadening 
and  the  shift  of  a  gold  Fermi  edge  as  a  function  of  spot  size,  laser  power,  and  emission  angle 
are  characterized  for  pulse  lengths  of  6  ps  and  6  eV  photon  energy.  The  results  are  used  as  a 
benchmark  for  an  A-body  numerical  simulation  and  are  compared  to  different  regimes  used  in 
photoemission  spectroscopy.  These  results  provide  an  important  reference  for  the  design  of  time- 
and  angle- resolved  photoemission  spectroscopy  setups  and  next- generation  light  sources. 
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I.  INTRODUCTION 


Electrons  traveling  in  vacnnm  to  a  measnrement  apparatns  will  repel  each  other  becanse  of 
their  mntnal  Conlomb  interaction.  If  a  clond  of  electrons  is  snfhciently  dense,  this  vacnnm 
space  charge  effect  can  alter  the  time,  energy,  and  spatial  spread  of  the  clond,  limiting 
experimental  resolntion  and  resnlting  in  systematic  measurement  errors.^®^ 

The  space  charge  effect  has  been  addressed  mostly  for  the  design  of  electron  guns  for 
time-resolved  electron  diffraction  experimentd^^  and  has  been  discussed  in  terms  of  N- 
body  numerical  simulations,  mean-held  models,  and  huid  dynamics.  Recently  the  space 
charge  effect  has  also  been  discussed  in  angle-resolved  photoemission  spectroscopy  (ARPES) 
experiments^,  where  it  has  been  shown  to  depend  primarily  on  four  parameters:  photon 
energy,  pulse  duration,  incident  beam  spot  size,  and  the  number  of  electrons  photo-emitted 
per  pulse.  These  previous  studies  explored  two  different  regimes:  i)  the  relatively  high 
photon  energy  (>  35  eV)  and  long  pulse  duration  (?«  60  ps)l^  regime;  and  ii)  the  very  low 
photon  energy  (two-photon  photoemission  with  3  eV  photons)  and  short  pulse  duration 
40  fs)  regime.®  The  shape  of  the  electron  cloud  is  very  different  in  these  two  regimes. 
In  the  hrst  case,  the  wide  distribution  of  electron  kinetic  energy  and  long  pulse  duration 
produce  a  cloud  extended  mainly  along  the  direction  normal  to  the  sample  surface,  while 
in  the  second  case  the  narrow  electron  kinetic  energy  distribution  and  short  pulse  duration 
produce  a  disk-shaped  electron  cloud  parallel  to  the  sample  surface.  The  former  case  of  a 
more  three-dimensional  cloud  requires  numerical  simulations  while  a  disk-shaped  cloud  can 
be  characterized  analytically  using  a  simple  power  law.®® 

With  the  recent  advent  of  laser-based  photoemission  experiments,  the  need  to  extend 
these  studies  to  an  intermediate  regime  becomes  imperative.  A  recent  study  of  the  space 
charge  effect  using  laser  based  ARPES  with  7  eV  photon  energy  and  10  ps  long  pulse  suggests 
that  the  space  charge  effect  plays  a  minor  role  in  this  regime  thanks  to  the  mirror  charge 
effect  that  nearly  cancels  the  space  charge  effect  in  the  1  to  10  ps  pulse  duration.l^ 

Here  we  expand  these  results  by  reporting  a  systematic  study  of  the  space  charge  effect 
on  polycrystalline  gold  in  laser-based  photoemission  spectroscopy  in  the  intermediate  regime 
of  moderate  pulse  lengths  (~  6  ps)  and  low  photon  energy  (6  eV).  We  hnd  in  the  limit  of 
small  spot  size  (~  20  /rm  diameter)  a  shift  of  the  Fermi  edge  on  the  order  of  2  meV  and  a 
broadening  larger  than  10  meV  starting  already  at  100  electrons  per  pulse.  More  generally. 
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we  find  that  a  single  power  law  cannot  describe  the  broadening  or  the  Fermi  edge  shift  as 
a  function  of  electrons  per  pulse  and  spot  size  in  the  various  ARPES  regimes.  To  explore 
the  non-trivial  relation  between  the  various  experimental  parameters,  we  modeled  our  data 
with  a  molecular-dynamics  model  which  employs  full  A^-body  numerical  simulations  of  the 
electron  propagation  from  the  surface  to  the  detector.®  This  model  is  able  to  reproduce 
quantitatively  the  space  charge  effect  seen  in  the  long  and  short  pulse  regimes  as  well  as  in 
the  intermediate  regime  explored  here.  The  large  range  of  applicability  of  this  model  allows 
one  to  identify  the  functional  dependence  of  the  spectral  shift  and  broadening  on  the  number 
of  electrons  per  pulse  and  spot  size.  We  found  that  one  of  the  key  parameters  affecting  the 
different  power  laws  observed  in  different  ARPES  regimes  is  the  maximum  kinetic  energy  of 
the  initial  photo-electron  distribution.  We  believe  that  this  study  will  be  crucial  for  guiding 
next  generation  high  resolution  setups. 

II.  EXPERIMENTAL  TECHNIQUE 

Experiments  were  performed  on  a  polycrystalline  gold  sample  scratched  in  air  and  inserted 
through  a  load  lock  into  a  UHV  analysis  chamber  at  low  temperature  (15-20  K).  Pressure 
remained  below  5  x  10“^^  Torr  throughout  the  experiment.  In  Fig.  1  we  show  a  diagram  of 
the  experimental  setup.  The  sample  geometry  is  detailed  in  the  inset.  The  spot  size  was 
characterized  by  using  a  wire  mesh  mounted  on  the  same  sample  holder  as  the  gold  sample. 

The  laser  system  is  a  cavity-dumped,  mode- locked  Ti: Sapphire  oscillator  (Coherent  Mira), 
pumped  with  a  6  W  frequency-doubled  Nd:YV04  laser.  The  oscillator  generates  ~  150  fs 
pulses  at  840  nm  at  a  repetition  rate  of  54.3/n  MHz  (n  =  1,  2,  3, . . . )  with  pulse  energies  as 
high  as  100  nJ.  The  photon  energy  of  the  laser  beam  was  frequency-quadrupled  via  cascaded, 
type-I  phase-matched  second  harmonic  generation  in  two  beta  barium  borate  (BBO)  crystals 
of  2  mm  and  5  mm  thickness.  The  resulting  210  nm  (5.9  eV)  pulse  was  estimated  to  be  ~  6 
ps  long  due  mainly  to  group  velocity  mismatch.  The  UHV  chamber  window  is  made  of  UV- 
grade  fused  silica.  We  adjusted  the  pulse  energy  by  tuning  the  power  applied  to  the  cavity 
dumper  Bragg  cell.  The  focusing  lens  (L3  in  Fig.  1)  has  a  300  mm  focal  length.  We  varied 
the  spot  size  by  defocussing  L3. 

Assuming  a  work  function  of  about  4.4  eV,  electrons  at  the  Fermi  energy  exit  the  sample 
with  1.5  eV  kinetic  energy,  corresponding  to  a  velocity  of  0.7  pm/ps.  With  a  pulse  lasting 
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FIG.  1:  (Color  online)  Laser-ARPES  experimental  setup.  The  focal  lengths  are:  Ll=10cm, 
L2=20cm  and  L3=30cm. 

6  ps,  the  initial  electron  spread  normal  to  the  surface  is  thus  about  4  pm.  With  a  spot 
size  diameter  of  about  20-100  pm,  the  initial  electron  cloud  is  a  flat  cylinder.  This  is  an 
intermediate  case  between  the  fs  laser-based  data  and  the  synchrotron-based  data.® 

The  integrated  photoemission  signal  was  measured  using  a  Phoibos  150  mm  hemispherical 
analyzer  from  Specs  with  a  0.2  mm  entrance  slit  and  4  eV  pass  energy,  yielding  an  energy 
resolution  of  4  meV.  The  sample  was  electrically  grounded  through  a  picoammeter,  which 
allowed  the  measurement  of  the  total  photo-electron  yield.  All  the  spectra  were  averaged 
over  18°  along  the  analyzer  slit  except  where  explicitly  stated  otherwise.  The  angle 
between  the  analyzer  and  the  laser  beam  is  hxed  at  55°  and  the  light  incident  on  the  sample 
is  s-polarized. 

Fig.  2  shows  the  characterization  of  the  pulse  energy  and  spot  size.  The  total  photo¬ 
electron  yield  as  a  function  of  the  laser  power  is  shown  in  Fig.  2(a).  The  repetition  rate 
was  set  to  5.4  MHz.  The  solid  line  is  a  linear  £t  giving  a  yield  of  1  x  10“'^  e“ /photon.  For 
comparison,  the  typical  yield  for  a  photon  energy  of  35  eV  measured  at  BL  10.0.1  at  the 
Advanced  Light  Source  (ALS)  is  about  0.2  e“/photon  (1  x  10^^  photons/s  with  300  nA 
yield)  The  difference  can  be  mainly  attributed  to  the  range  of  electronic  states  probed:  6 
eV  photons  only  excite  states  in  the  An  5d-6s  valence  band  with  binding  energy  less  than 
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FIG.  2:  (Color  online)  Spot  size  measurements,  (a)  Photo-electron  yield  as  a  function  of  laser 
power,  (b)  Typical  x  and  y  scan  of  the  integrated  photoemission  intensity  from  the  W  mesh.  The 
y  scan  is  vertically  shifted 

about  1.5  eV,  whereas  35  eV  photons  can  excite  the  entire  valence  band.^  Additionally,  in 
the  latter  case  secondary  electrons  amplify  the  photo-electron  yield  noticeably.^ 

The  spot  size  was  measured  by  scanning  a  100x100  lines/inch  tungsten  mesh  at  the 
sample  position  while  measuring  the  integrated  photoemission  intensity.  The  mesh  wire 
diameter  was  25  pm.  The  6  angle  (see  Fig.  1)  was  set  to  55°  during  the  mesh  scan  so  that 
the  X  and  y  scanning  directions  were  perpendicular  to  the  laser  beam.  Fig.  2(b)  shows  a 
typical  scan  along  the  x  and  y  direction  obtained  from  this  procedure.  We  then  estimated 
the  spot  size  by  deconvolving  the  peak’s  FWHM  using  the  known  mesh  wire  diameter.  The 
reported  diameter  is  the  square  root  of  the  x  and  y  FWHM  product,  corresponding  to  the 
diameter  of  a  disc  of  equal  area.  The  ratio  between  the  FWHM  along  x  and  y  is  not  constant 
but  was  smaller  than  3.  We  note  that  with  an  incident  Gaussian  beam  of  2  mm  diameter, 
the  diffraction  limited  spot  size  with  a  300  mm  focal  length  lens  is  about  15  pm.  For  the 
FWHM  along  the  x  direction,  we  used  a  corrected  wire  diameter  to  account  for  the  fact  that 
only  the  part  of  the  wire  facing  the  analyzer  could  contribute  to  the  measured  signal.  It 
can  be  shown  that  the  effective  wire  diameter  in  this  case  is  25  x  (1  -f  cos(55°))/2  =  20  pm. 
The  small  difference  in  peak  separation  between  the  x  and  y  scan  is  due  to  mesh  spacing 
inhomogeneity. 

Fig.  3  shows  typical  integrated  photoemission  spectra  from  the  gold  sample  at  various 
laser  powers  showing  the  Fermi  level  broadening  and  shift.  The  spot  size  diameter  was  18 
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FIG.  3:  (Color  online)  Fermi  edge  of  the  An  sample  as  a  function  of  the  number  of  photo-electrons 
per  pulse  (A^).  The  spectra  are  normalized  and  shifted  vertically.  The  solid  lines  are  fits. 

/rm.  The  edges  of  these  spectra  have  been  £t  to  a  function  of  the  form: 

f{E)  =  a  +  b-  mm{E  -  Eg,  0)  +  c  •  erfc  •  (1) 

The  parameters  a,  b,  c,  Eq,  and  FWHM  are  htting  parameters  and  E  is  the  electron  kinetic 
energy.  The  hrst  two  terms  account  for  the  background,  which  includes  an  offset  and  a  linear 
contribution  for  E  <  Eq.  The  third  term  is  the  complementary  error  function  erfc(x), 

erfc(x)  =  1  —  erf(x)  =  1 - ^  (2) 

VTT  Jo 

which  is  commonly  used  to  fit  Fermi-Dirac  distributions  convolved  with  experimental  reso¬ 
lution.  The  characteristic  FWHM  is  defined  to  be  the  full-width-half-maximum  value  of  the 
associated  normalized  Gaussian. 

The  top  of  the  panel  shows  a  sharp  Fermi  edge  at  low  power  with  a  FWHM  of  7.3  meV, 
consistent  with  the  experimental  energy  resolution  (4  meV  from  the  analyzer  and  4-5  meV 
from  the  laser  bandwidth)  and  the  sample  temperature  (about  15  K).  The  total  photo¬ 
electron  yield  measurement  indicated  an  average  of  7  electrons  per  pulse.  As  the  laser  power 
is  increased,  a  greater  number  of  photo-electrons  (N)  are  emitted  per  pulse.  This  results 
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in  a  stronger  space  charge  effect,  and  the  edge  becomes  broader  and  shifts  to  higher  kinetic 
energy. 

The  shift  to  higher  kinetic  energy  is  due  to  the  fact  that  the  electrons  at  the  Fermi 
energy  have  the  highest  initial  kinetic  energy,  and  therefore  the  Coulomb  repulsion  of  slower 
photo-electrons  can  only  push  them  towards  higher  kinetic  energy  (neglecting  any  mirror 
charge  effects).  We  note  that  within  the  laser  power  range  used  here,  the  erfc(x)  fit  function 
captures  the  experimental  line  shape  very  well.  This  was  not  the  case  at  higher  laser  power 
(not  shown). 


III.  EXPERIMENTAL  RESULTS 

In  Fig.  4  we  report  the  shift  of  the  Fermi  edge  Aef,  and  its  broadening  F,  as  a  function  of 
number  of  electrons  per  pulse  (iV)  for  several  values  of  the  spot  size  diameter  {D).  Here  Aef 
is  defined  by  the  difference  between  the  measured  edge  position  (Eq)  and  the  Fermi  energy 
{Ef)  while  F  is  defined  as  the  edge  FWHM  deconvolved  from  the  sample  temperature  and 
experimental  resolution.  The  F  error  bars  are  mainly  estimated  from  fit  standard  deviation 
and  temperature  measurement  uncertainty  (estimated  to  be  about  3  K).  The  Aef  error  bars 
are  the  standard  deviations  plus  0.5  meV.  The  error  on  N  was  estimated  to  be  ±(30-|-10%). 

We  observe  that  F  and  Aef  increase  with  increasing  N  and  that,  as  expected,  this  effect 
becomes  more  severe  as  the  spot  size  decreases.  The  data  from  Ref.  [H]  taken  in  a  different 
regime  {hu  =  35  eV  and  pulse  length  =  60  ps),  corresponding  to  the  largest  D  value  shown 
in  the  figure  (see  diamond  markers),  fall  in  line  with  this  trend.  To  determine  the  relation 
between  F,  Aef,  and  N  we  fit  the  data  with  the  smallest  D  value  using  the  same  power 
law  as  the  one  reported  in  previous  studies.®  It  is  clear  that  neither  the  linear  (see  dashed 
line)  nor  the  y/N  dependence  (dotted  line)  can  describe  effectively  the  data  here  reported 
in  the  ps  regime.  To  model  the  transition  from  the  linear  to  the  y/N  regime  we  fit  the  data 
of  Fig.  4  with  generalized  exponents: 


T{N)  =a^x 

(3) 

EFiN)  =  X 

(4) 

In  Figs.  5(a)  and  (b),  we  report  the  dependence  of  the  fit  coefficients  a  and  j3  on  D, 
respectively.  Remarkably,  (3  is  essentially  independent  of  D  for  both  F  and  Aef-  More 
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FIG.  4:  (Color  online)  Space  charge  effect  as  a  function  of  the  number  of  electrons  per  pulse  (A^) 
for  various  spot  diameters  {D).  (a)  Broadening  (F)  of  the  Fermi  edge  as  a  function  of  D  and  N. 
The  dashed  line  shows  a  linear  fit  and  the  dotted  line  a  '/N  fit  to  the  triangular  markers,  (b)  Shift 
of  the  Fermi  edge  {Aef)  as  a  function  of  N  for  various  values  of  D.  The  dashed  lines  show  a  linear 
fit  of  the  triangular  markers.  The  solid  lines  are  fits  to  a  power  law.  The  data  set  with  diamond 
shaped  markers  was  obtained  at  the  ALS  with  35eV  photon  energy  and  was  adapted  from  Ref.  [9l 
The  y-axis  error  bars  include  the  errors  on  the  fit  and  the  instrumental  resolutions.  The  x-axis 
error  bar  is  estimated  from  the  sample  current  measurement  resolution.  The  spectra  are  averaged 
over  ~18°  along  the  analyzer  slit  with  9  =  55°. 

8 


FIG.  5:  (Color  online)  Plot  of  the  the  fitting  coefficients  a  and  /3  as  a  function  of  D.  The  solid 
lines  in  panel  (a)  show  the  average  value  of  /3r  and  Paef-  The  solid  and  dashed  lines  in  panel  (b) 
are  fits  to  a  power  law. 


specifically,  for  the  broadening  we  found  =  0.7,  which  is  an  intermediate  value  between 
the  value  (1)  reported  for  synchrotron  pulsed  and  the  value  (0.5)  reported  for  fs  pulses.® 
For  the  Fermi  edge  shift  we  found  Paef  =  1-1.  In  both  cases  the  prefactor  a  increases  with 
decreasing  D.  To  characterize  the  space  charge  effect  dependence  on  D  more  quantitatively, 
we  £t  a  with  a  function  of  the  form 


a{D) 


C 

1^  ’ 


(5) 


where  C  and  A  are  fitting  parameters.  In  a  log-log  plot,  a{D)  is  a  straight  line  with  slope 
A.  By  substituting  a{D)  into  Eq.  (|^,  we  hnd  that  the  space  charge  effect  dependence  on 
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FIG.  6:  (a)  F  and  as  a  function  of  the  emission  angle  cj).  (Aef)  =  77  meV  and  (F)  =  79 

meV  are  estimated  from  the  fit  of  the  angle-integrated  spectrum,  (b)  Photoemission  intensity  as  a 
function  of  cj).  The  solid  lines  are  fits  with  the  function  /i((/))  =  const  x  cos^ (</>). 


the  number  of  photo-electrons  N  and  spot  diameter  D  is  of  the  from: 

f(N,D)  =  Cx^. 


More  specihcally  we  found  for  the  Fermi  edge  shift: 

ATi-i 

^EFiN,D)  =  0.5  X  , 


(6) 

(7) 


and  for  the  Fermi  edge  broadening: 

NO.7 

r(/v,r>)  =  i3x  — ,  (8) 

where  as  in  all  following  formulas  Aef{N,  D)  and  F(iV,  D)  are  in  meV  and  D  is  in  pm. 

To  quantify  the  angular  dependence,  the  angle-resolved  results  are  shown  in  Fig.  6.  The 
emission  angle  (0)  is  along  the  same  axis  as  6  with  0  =  0°  corresponding  to  normal  emission. 
The  angular  acceptance  of  the  analyzer  was  ~  25°  along  the  entrance  slit.  Panel  (a)  shows  the 
strong  dependence  of  Aff  and  F  from  0.  We  show  in  panel  (b)  the  photoemission  intensity 
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measured  at  ~  200  meV  below  the  Fermi  energy.  As  a  guide  for  the  eye,  we  £t  Aef{4>),  r(0) 
and  the  photoemission  intensity  with  empirical  functions  of  the  form  /i(0)  =  const  x  cos^  0. 
An  offset  is  added  in  the  case  of  the  photoemission  intensity. 

The  strong  correlation  between  the  Fermi  edge  broadening  and  the  photoemission  inten¬ 
sity  suggests  that  the  angular  dependence  of  the  space  charge  effect  is  mainly  due  to  the 
photo-electron  emission  angular  distribution.  The  effect  of  the  angular  integration  over  the 
analyzer  acceptance  angle  is  different  for  Aef  and  F.  While  the  shift  of  the  angle-integrated 
spectrum  {{Aef)  =  77  meV)  is  close  to  the  average  value  of  AEF{(t>)  over  0,  the  broadening 
of  the  angle-integrated  spectrum  ((F)  =  79  meV)  is  significantly  larger  than  the  average 
value  of  F(0).  This  shows  that  the  broadening  due  to  the  angular  dependence  of  Aef{4>) 
can  be  larger  than  just  the  kinetic  broadening.  This  angular  distribution  of  the  space  charge 
effect  has  a  direct  effect  on  the  data  reported  here  since  the  spectra  were  integrated  over 
^  18°  with  0=55°. 

Since  we  reported  the  space  charge  effects  as  a  function  of  N,  measurements  at  different 
emission  angles  will  yield  different  values  for  the  coefficient  a.  This  should  result  in  a 
different  value  of  C.  The  exact  contribution  of  the  angle  integration  window  and  emission 
angle  on  the  exponents  0  and  A  is  more  difficult  to  estimate. 


IV.  SIMULATION  RESULTS 

The  numerical  simulations  are  based  on  a  representative  electron  energy  distribution  with 
a  sharp  Fermi  edge  at  1.56  eV.  The  pulse  duration  was  assumed  to  be  rectangular  with  a 
duration  of  6  ps.  The  simulations  were  repeated  at  least  10  times  for  each  parameters 
setting.  Aef  and  F  are  extracted  with  the  same  edge  fit  fnnction  with  the  exception  of  the 
parameters  a  and  b  being  set  to  0.  The  electron  emission  was  assnmed  to  follow  a  cos^(0) 
distribntion  following  the  experimental  results  of  Fig.  6  and  only  electrons  with  an  emission 
angle  between  46°  and  64° are  integrated  in  order  to  mimic  the  experimental  condition.  The 
nnmber  of  electrons  per  pulse  N  varied  between  200  and  4000.  Spot  diameters  between  18 
fim  and  359  /rm  were  used. 

Fig.  7  shows  the  simulated  broadening  F  and  shift  Aef  of  the  Fermi  edge  as  a  function 
of  the  nnmber  of  electrons  per  pulse  N  for  different  spot  diameters  D.  Fig.  8  shows  F  and 
Aef  a  function  of  the  spot  diameter  D  for  different  numbers  of  electrons  per  pulse  N. 
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FIG.  7:  (Color  online)  Simulation  of  the  space  charge  effect,  (a)  F  and  (b)  IS.ef  as  a  function 
of  N  for  various  D  value.  The  error  bars  result  from  multiple  simulation  runs  and  represent  the 
standard  deviations.  The  solid  lines  are  fits  to  a  power  law. 

We  fit  the  results  for  F  and  Aef  with  the  same  power  law  as  in  the  experimental  case  and 
show  the  comparative  results  in  Table  1.  The  exponent  values  are  the  average  value  of  all 
the  htting  results  (not  shown). 

The  broadening  and  the  shift  of  the  Fermi  edge  show  a  good  overall  agreement  between 
simnlation  and  experiment.  Small  systematic  deviations  are  observed  in  the  small  effect 
regime.  In  particnlar  for  N  <  600  and/or  D  >  200  /im. 
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FIG.  8:  (Color  online)  Simulation  of  the  space  charge  effect,  (a)  F  and  (b)  as  a  function 

of  D  for  various  N  values.  The  error  bars  result  from  multiple  simulation  runs  and  represent  the 


standard  deviations.  The  solid  lines  are  fits  to  a  power  law. 


V.  DISCUSSION 


In  the  case  of  the  disk-shaped  photo-electron  cloud  produced  by  a  fs  pulse,  the  broadening 
due  to  the  space  charge  effect  was  proposed  to  be  of  the  forrrP 


r  oc 


(9) 


On  the  other  hand,  for  the  60  ps  synchrotron  pulse,  the  broadening  follows  a  linear  relation 
with  N .  The  F  oc  N^''^  for  the  simulation)  dependence  we  report  for  this  intermediate 
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Simulation 

Experiment 

Aef{N,D)  = 

/YO.79 

3-lX  ^0.85 

0'5x  jji  i 

T{N,D)  = 

jyO.83 

5-5x  £)i.o2 

NO.7 

13  X  '^1.2 

TABLE  I:  Summary  of  the  measured  and  simulated  exponent  values. 

regime  of  laser-based  ARPES  with  6  ps  pulses  is  between  these  two  values.  However,  these 
three  regimes  cannot  be  solely  compared  on  a  pulse  duration  basis.  The  fs  pulse  analysis 
is  based  on  a  surface  state  peak  (86  meV  FWHM)  measured  at  the  band  minimum  in  the 
Brillouin  zone  center.  In  this  case,  the  majority  of  the  photo-electrons  have  a  higher  kinetic 
energy.  This  will  affect  the  space  charge  effect  in  two  ways:  i)  the  peak  position  will  be 
shifted  to  lower  kinetic  energy;  ii)  the  momentum  resolution  will  decrease  resulting  in  a 
broader  spectrum  shifted  towards  higher  kinetic  energy.  While  the  total  shift  in  energy  ends 
up  being  negligible,  the  contribution  of  the  momentum  broadening  on  T  is  not.®  This  makes 
a  direct  comparison  difficult  but  further  comparative  analysis  can  be  done  with  our  model. 

According  to  the  simulation  results,  the  nonlinearity  T  oc  A^^^  with  /dr  7^  1  seems  to 
depend  largely  on  the  initial  photo-electron  distribution.  More  specifically,  the  simulations 
show  that  the  exponent  /dr  decreases  as  the  maximum  electron  kinetic  energy  decreases. 
This  exponent  is  around  1  for  high  kinetic  energies,  decreases  to  /dr  =  0.7  —  0.8  in  the 
intermediate  regime  presented  here  with  a  maximum  electron  kinetic  energy  of  about  1.5 
eV,  and  decreases  down  to  /dr  =  0.5  for  the  fs  pulse  experiment  where  the  maximum  electron 
kinetic  energy  was  around  1.1  eV.®  We  found  also  that  the  simulation  shows  no  indications 
that  the  pulse  durations  affect  this  behavior  of  the  exponent  /dr. 

Similarly,  we  hnd  the  broadening  dependence  on  spot  size  to  be  proportional  to  D~^'^ 
{D~^  for  the  simulation),  compared  to  D~^'^  the  fs  pulse  data,  though  more  experimental 
data  are  required  here  to  validate  the  trend. 

It  is  interesting  to  note  also  that  the  ratio  of  the  two  coefficients  /dr  and  Ar  is  not  1  as 
for  the  fs  pulse  case,  but  closer  to  0.6-0. 8.  Since  the  error  bars  on  A  are  relatively  large, 
fits  with  A  =  /d  are  shown  with  a  dashed  line  in  Fig.  5  (b).  From  this  ht  it  is  clear  that 
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additional  data  points  are  needed  to  conclude  on  the  correct  value  of  the  ratio  between  A 
and  13.  Nevertheless,  these  results  set  already  a  limit  in  the  range  in  which  the  analytical 
model  proposed  in  Ref.  El  is  valid. 


VI.  CONCLUSIONS 

In  conclusion,  we  report  the  measurement  of  the  space  charge  effect  in  the  few-ps  laser 
pulse  regime.  We  observe  a  broadening  of  the  spectral  feature  exceeding  10  meV  and  an 
energy  shift  exceeding  2  meV  with  pulses  of  as  few  as  a  100  photo-electrons  for  small  spot 
sizes  (~  20  pm  diameter).  For  larger  spot  sizes  (~  165  pm  diameter)  similar  effects  are 
observed  only  as  the  number  of  photo-electrons  per  pulse  reaches  600  to  800.  We  found 
that  neither  the  linear  model  that  hts  synchrotron  data  nor  an  \/N  model  that  hts  fs  pulses 
describe  our  data  effectively,  and  we  introduced  a  simple  phenomenological  model  with 
generalized  exponents  that  characterizes  the  broadening  and  energy  shift  of  the  spectra  as 
a  function  of  spot  size  and  number  of  photo-electrons  over  the  entire  range  probed.  The 
results  were  used  as  a  benchmark  for  an  iV-body  numerical  simulation  with  a  good  overall 
agreement.  In  turn  the  simulation  allowed  us  to  explain  the  difference  observed  in  different 
ARPES  regimes  in  terms  of  initial  electron  kinetic  energy.  These  results  represent  a  reference 
for  the  design  and  interpretation  of  photoemission  experiments  with  next-generation  light 
sources,  such  as  free-electron  lasers  and  table-top  laser  sources. 
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